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INTRODUCTION 

Problems  Involving  network  analysis  arise  In  many  Industrial  contexts.  For  example, 
interest  often  attaches  to  methods  for  characterizing  the  behavior  of  firewater  safety  systems 
in  refinery  operations.  Flows  and  resulting  pressures  need  to  be  ascertained  under  various 
patterns  of  input  and  withdrawal.  Under  certain  conditions  standard  approaches  such  as  Hardy 
Cross*  methods  and  extensions  of  the  usual  Gaiisa-fieW^I  techniques,  have  been  found  difficult 
or  impossible  to  apply  in  analyzing  such  systems.  The  nonlinear  character  of  the  branch  head 
losses  may  cause  lack  of  convergence  of  the  solutions  or  require  recourse  to  supplementary 
analyses  of  a  complicated  character. 

•  Certain  kinds  of  network  problems  have  proved  amenable  tc  the  usual  |ormuxiffo'ts  pf 
linear  programming,^  flrne  it  is  -preyosed-t^-t^iiiTg  a  wanewhat  different  appfdflCh^  order  to 
provide  a  more  general  method  4  attack  and  to  tin  we,  moreover;  -in  •  W*y  whi#»  provides 
access  to  the  methods  4eewibefHtrf4i  for  treeing  nonlinear  problems  as  well.  The  essential 
ideas  are  as  follows:  (1)  re  statement  of  the  original  problem  in  terms  of  minimizing  a  sepa¬ 
rable  convex^ functional  subject  to  a  system  of  linear  equations  expressing  the  conservation  of 
current  entering  and  leaving  nodes  by  branches;*  (2)  u*Uizati,on  ct  piecewise  linear  functionals 
to  secure  suitable  approximations  to  the  convex  functiona}s;|(3)  using  the  "bounded  variables 
tech;dques**^cd  linear  programming  in  order  to  avoid  expanding  the  number  of  vectors  In  an 
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^See  M 

2See,  e.g.,  W  Prager,  [9]  and  the  references  cited  therein. 

3 By  "separable  convex"  function  is  meant  one  which  is  the  sum  of  convex  functions  of  one 
variable  each 

4i.e.,  it  is  assumed  that  leakages,  if  any,  may  be  ignored. 
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•clive  basis  snd  aiso  to  avoid  expanding  the  number  of  dl/ferent  coefficient  vector*  which 
require  explicit  expression.  The  resulting  problem  is  thereby  converted  to  an  equivalent  one 
which  may  be  called  the  polygonal  approximation  to  convex  programming  over  the  'Incidence 
matrix" of  a  connected  network.  Special  properties  are  then  secured  which  simplify  analysis 
and  representation,  and  provide  a  basis  for  systematic  elaboration, 

INCIDENCE  MATSZ^ 

In  engineering  analyses  of  networks.  It  is  customary  to  designate  one  of  the  two  pos¬ 
sible  directions  of  flow  in  a  branch  as  positive.  The  currert  is  then  allowed  to  have  either 
positive  or  negative  values.  Positivity  and  negativity  are  then  respectively  associated  with 
flow  in  the  positive  direction  or  its  opposite. 

Topological  analysis  utilizes  instead  the  closely  related  notion  of  "incidence"  of 
branches  (or  links)  on  nodes.  Thus,  the  convention  may  be  -  ■footed,  for  example,  that  a  branch 
entering  a  node,  is  positively  incident  on  that  node  a  branch  leaving  a  node  is  negatively  inci¬ 
dent  on  it,  and  a  branch  which  does  not  meet  a  node  is  not  incident  on  it.  The  incidences  of 
branches  on  nodes  may  then  be  described  by  an  incidence  matrix  consisting  of  numbers  ( 
representing  the  incidence  of  the  jth  branch  on  the  1th  node  as  follows: 


(1) 


(  1  |  if  branch  J  is  positively  incident  on  node  i . 
<-l  >  "  "  J  "  negatively  ”  "  "  i. 

0  "  "  j  "  not  "  "  "  l , 


Since  each  branch  is  incident  on  precisely  two  nodes  each  column  of  thss  matrix  will 
contain  precisely  two  nonzero  entries— a  plus  one  and  a  minus  one.  If  designates  the  cur¬ 
rent  in  the  branch  the  net  Influx  of  current  Into  the  i1*1  node  is  evidently 


(2) 


\  £ij  qj  ’ 


where  J  ranges  over  all  branches  of  the  network.  If  the  efflux  out  of  the  network  at  this  node 
Is  Ej  then  this  conservation  of  current  at  nodes  (Kirchhoff's  node  law)  may  be  written  as 


(3) 


^  flJ  qj  =  El’ 


1=1,2, 


Thus,  to  express  the  conservation  conditions,  It  Is  necessary  only  to  form  the  Incidence  matrix 
column  by  column. 

BASES  OF  COLUMN  VECTORS  FOR  AN  INCIDENCE  MATRIX 

As  stated  In  (1)  above,  each  column  vector  of  such  a  matrix  contains  a  plus  one  and  a 
minus  one  and  is  uniquely  associated  with  a  branch  of  the  network.  The  plus  one  corresponds 
to  the  node  at  the  head  of  the  branch  and  the  minus  one  corresponds  to  the  node  at  the  tall  of 
the  branch.  A  basis  of  column  vectors  of  the  matrix  Is  a  collection  of  column  vectors  (or 
branches)  such  that,  (1)  every  other  column  vector  can  be  expressed  as  a  linear  combination 
of  this  collection,  and  (2)  the  number  of  vectors  In  this  collection  is  the  least  with  which  this 
task  can  be  accomplished.  Clearly,  a  basts  must  lncluae  at  least  one  vector  incident  (posi¬ 
tively  or  negatively)  on  each  node  of  the  system.  Otherwise,  It  would  not  be  possible  to  expiess 
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(in  terms  of  the  basis)  a  v«  t«>r  associated  with  a  branch  Incident  on  ttut  node,  To  express  a 
column  vector  (branch  or  node  pair)  In  terms  of  the  basis.  It  1*  necessary  only  to  start  with 
the  head  node  of  this  branch  and  trace  a  path  alor^  branches  associated  with  the  basis  to  the 
tall  node  of  this  branch.  The  corresponding  algebraic  process  consists  of  adding  or  subtract¬ 
ing  the  basis  vectors  associated  with  the  branches  thus  traversed. 

Consider,  for  example,  the  process  of  moving  from  one  branch  to  another  across  a 
node.  Both  branches  arc  Incident  on  this  node.  If  both  have  the  same  Incidence,  subtracting 
them  will  give  a  column  vector  with  precisely  one  "plus  one"  and  one  "minus  one”  associated 
with  the  two  end  (or  exterior)  nodes  of  these  two  branches;  11  the  branches  are  oppositely 
Incident  on  a  node,  an  addition  will  be  performed  to  attain  a  vector  corresponding  to  tie  end 
nodes.  Thus,  traversal  of  the  path  from  node  to  node  is  associated  with  addition  and  subtrac¬ 
tion  of  vectors  so  that  as  any  particular  node  in  the  traversal  Is  reached  the  cumulant  vector 
has  nonzero  entries  (of  opposite  sign  and  unit  magnitude)  corresponding  precisely  to  the  ei  d 
nodes  of  the  completed  portion  of  the  path. 

The  collection  of  branches  corresponding  to  the  basis  can  contain  no  closed  path  (or 
cycle).  Otherwise  It  would  be  possible  to  express  one  of  the  branches  In  the  cycle  in  terms  of 
the  others  by  the  process  just  described.  Every  expression  containing  the  vector  of  this 
branch  could  then  be  replaced  by  an  expression  not  containing  It.  The  supposed  basis  would 
then  contain  more  than  the  minimum  number  of  vectors  required  and  would  therefore  not  be  a 
basis. 

A  connected  system  of  branches  and  nodes  containing  no  cycle  is,  In  topology,  called  a 
tree."  Thus,  a  basis  for  the  incidence  matrix  of  a  connected  network  consists  of  the  column 
vectors  associated  with  a  tree  which  contains  every  node  of  the  network.  An  easy  method  is 
thereby  provided  both  for  securing  a  basis  and  fo?'  expressing  all  vectors  in  terms  of  It. 

DETERMINING  THE  "EVALUATORS" 

The  "evaluators"— the  w  vector6 *— may  be  determined  In  a  manner  Blmilar  to  the  way 
In  which  the  "row-column  numbers"  are  determined  in  distribution  type  models  of  linear 
programming. 1  In  the  model  which  Is  belrg  developed  Is  associated  with  the  1th  equation. 
The  value  of  any  one  such  Wj  may  be  selected  arbitrarily  because  the  system  of  equations  (3) 
has  the  same  degree  of  linear  dependence  as  (he  standard  distribution  model.8  The  other  u>. 
are  then  quickly  determined  since  each  basis  branch  consists  of  a  node  pair  and  the  calculations 
are  initiated  at  node  1  and  proceed  away  from  It  to  all  the  other  nodes  of  the  network  along  the 
appropriate  basis  branches  (hence  column  vectors). 

EXTREMAL  PRINCIPLE  AND  THE  FUNCTIONAL 

In  the  case  of  flow  problems  to  be  considered  here,  It  Is  assumed  that  to  each  branch  J 
there  corresponds  a  resistance  function  r|  (qj)  such  that 


6Any  solution  of  Wl  B  *  c^,  where  B  is  the  basis  at  the  current  iteration.  3ee  [7],  for 
example. 

'In  the  following  discussion  some  familiarity  with  linear  programming  is  assumed.  See, 
e.g.,  [3]  and  [5]  for  elaboration. 

8This  may  be  seen  from  the  fact  that  the  number  of  branches  in  a  basis  tree  is  one  less 
than  the  number  of  nodes  In  physical  terms  this  corresponds  to  the  fact  that  the  voltage  dif¬ 
ferences  rather  than  absolute  voltages  are  determined  by  the  currents  in  a  network. 
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r.  fql  dq 


is  a  convex  function  of 


(5) 


For  example,  In  the  firewater  system  problem 


rj  (q) 

and 


spn  (q)  x 


Rj(q) 


J 

k+  1 


k* 1 


k  '  1 


The  convex  programming  problem,  or  extremal  principle,  characterizing  the  network 

flow  is 

min.  R,  (q4) 

J  J  ] 

(6)  subject  to 

^  ( 1J  =  ^1  ’  1  =  1 .  2,  . . . 

It  may  be  established9  that  any  optimal  solution  to  this  problem  automatically  satisfies  the 
"head-loss,"  or  Kirchhoff  cycle  conditions— viz,  the  algebraic  sum  of  the  head-losses  around 
any  closed  loop  consisting  entirely  of  current  carrying  branches  is  zero.  This  principle  thus 
achieves  an  essential  simplification  of  the  problem  because  in  arriving  at  a  solution  it  is  only 
necessary  to  consider  individual  branches  with  the  highly  nonlinear  loop  conditions  being 
transferred  to  the  functional10  in  a  much  less  recondite  form. 


APPROXIMATION  OF  THE  FUNCTIONAL 

To  any  desired  degree  of  approximation  each  Rj  (q)  may  be  replaced  by  a  piecewise 
linear  function  so  that  the  total,  £  R,  (q^),  may  also  be  approximated.11  The  resulting  problem 

may  then  be  replaced  by  an  equivalent  linear  programming  problem  using  the  bounded  variables 
techniques  in  order  (a)  to  avoid  proliferation  of  the  number  of  vectors  in  the  active  basis  and 
(b)  to  avoid  increase  in  the  number  of  structurally  distinct  vectors.12 


9The  authors  wish  to  thank  the  referee  for  stimulating  i  sufficiently  general  proof  of  this 
(which  will  be  published  elsewhere)  based  on  extending  the  Kuhn-Tucker  theorem  for  this  situa¬ 
tion.  He  pointed  out  that  Duffin,  Birkhoff,  and  Diaz  required  rj  (q)  continuous  and  strictly 
increasing  in  their  proofs. 

10lt  should  be  noted  that  the  objective  he  re  is  to  obtain  a  simulation  and  that  this  is  accom¬ 
pli  shed  by  optimization  with  re  suiting  gains  in  efficiency  for  the  desired  simulation.  For  mathe¬ 
matical  proofs  of  this  principle  under  heavier  restrictions  than  are  necessary.  seeflOland 

[H], 

1  Additional  research  which  eliminates  the  need  forthis  approximation  stepwdll  be  reported 
in  a  subsequent  paper.  The  method  to  be  r epo rted  con sist s  essentially  of  implicit  refinement  to 
the  limit  of  the  approximation  process  and  of  short-hand  rules  so  that  no  approximation  needs 
to  be  dealt  with  explicitly,  This,  of  course,  replaces  the  linear  iterations  with  a  nonlinear 
counterpart. 

1 2See  [4]  for  a  general  description  of  this  method  and  others  for  approximat ing  a  nonlinear 
functional  by  a  piecewise  linear  one. 
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In  detail  i«< 


qj  *  V  ■  qj' :  V  ■  qj  0 


ar>d  further  let 

(7) 


r  q( 


k  J' 


'  A  ‘ 
k  >k 


0  5  A 


jk  ■*  \ 


0  *  '  V  • 


Then,  the  approximating  problem  to  (6)  assumes  the  form  12 

|Pjk+  Ajk  +  pjk  Ajk 


min.  12 

j,k  L 


subject  to 

(8) 


L 


'« [k(a'*k 


jk  -  Ajk  =  Ei 


0  <  V*  *  V 


0  £  V  *  V 


vrhere  i,  j,  and  k  are  each  assigned  their  respective  ranges. 


METHODS  OF  COMPUTATION  AND  ILLUSTRATION 

The  method  of  computation  on  the  network  graph  may  be  summarized  as  follows  (the 
analogous  steps  on  a  tableau  will  be  clear  from  the  illustration).  One  starts  with  a  solution  of 
the  node  conditions  admitting  only  branches  which  form  a  basis.  The  w -vector  Is  then  calcu¬ 
lated  using  the  active  basis,  and  from  it,  the  Zj  is  simply  u  ^  -  u.'t ,  where  node  h  is  the 
"head"  of  branch  j  and  node  t  is  its  "tall"  as  determined  by  its  orientation.  These  are  com¬ 
pared  as  usual  with  the  c^  (in  equation  (8)  the  Pj^)  to  determine  optimality  or  nonoptlmality. 

If  nonoptimal  the  vector  to  "come  in"  to  the  basis  is  expressed  in  terms  of  the  active  basis 
simply  by  going;  from  the  tail  to  the  head  of  the  come  In  branch  on  active  basis  branches. 

Those  branches  whose  orientation  is  opposite  to  this  direction  of  travel  will  have  their  "cur¬ 
rent"  increased  in  the  new  solution,  those  with  same  orientation  will  have  it  decreased.  In 
other  words,  the  special  structure  of  this  problem  is  used  to  generate  the  information  neces¬ 
sary  to  a  simplex  (or  dual)  procedure  at  each  stage  rather  than  using  the  general  algorithm 
and  storing  additional  information. 

The  example  of  Figure  1  may  be  used  to  provide  a  simple  Illustration  for  the  computa¬ 
tional  routine  to  be  employed.  In  this  Figure  an  influx  of  10  units  is  being  put  into  the  system 
at  node  1  and  effluxes  of  2,  5,  and  3  units,  respectively,  are  being  drawn  from  nodes  5,  7,and9. 
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Figure  1  -  Network  Graph 

The  coefficients  (here  constants)  for  the  variables  which  enter  Into  the  functional  are  shown 
Inside  the  diagram  opposite  each  arrowhead.  13  a  physical  analogy  for  these  coefficients  is 
found  In  the  resistances,  or  htad  losses,  over  these  branches  In  the  circuit  which  are  here 
assumed  to  be  constants.  The  resistances  are  also  assumed  to  be  the  same  Irrespective  of  the 
direction  of  flow. 14  Finally  the  branch  2,  5  Is  assumed  to  be  characterized  by  a  resistance  of 
one  up  to  4  units  of  current  and  a  resistance  of  two  for  greater  currents.  This  can  be  visual¬ 
ized  most  easily  by  inserting  two  branches  between  nodes  2  and  5,  designating  the  first  as  2^5, 
to  indicate  that  It  Is  limited  in  the  amount  of  flow  it  will  admit  and  the  second  by  2,  5  to  Indicate 
the  absence  of  such  an  upper  bound  on  the  possible  values  of  the  flows  through  this  (the  latter) 
branch.  The  objective  then  is  to  characterize  the  current  flows  by  minimizing  the  value  of  the 
functional  overall  branches  of  the  network,  taking  account  of  the  fact  that  2,1  is  limited  to  a 
capacity  of  four  units,  as  indicated  by  the  value  4  ]  alongside  Its  arr  owhead,  and  that  further 
flows  along  the  branch  will  encounter  the  increased  resistance  associated  with  2,  5. 

Table  1  contains  the  corresponding  incidence  matrix  with  nodes  entered  in  the  stub  and 
branches  at  the  head.  Each  column  Is  composed  of  the  values  as  described  In  section  2, 
above.  The  external  Influxes  and  effluxes  for  the  system  are  noted  at  the  right  of  the  matrix  In 
the  column  labelled  "stipulations."  Each  row  of  the  matrix  then  represents  an  equation  which 
Is  a  Kirchhoff  node  condition  to  describe  the  conservation  of  current.  At  the  top  of  each  branch 
Is  a  value,  c j ,  to  be  entered  Into  the  functional  and  drawn  from  the  amounts  entered  at  the 
arrowheads  on  the  branches  of  Figure  1. 

Computation  may  proceed  by  means  of  either  the  Table  or  the  Figure.  For  purposes  of 
exposition,  both  devices  will  be  used. 

First  an  active  basis  Is  selected.  To  guarantee  that  a  basis  is,  In  fact,  selected  It  Is 
necessary,  as  previously  described,  only  to  choose  branches  in  such  a  manner  that  all  nodes 
in  the  system  are  attained  without  loops.  Hence  computations  may  commence  by  selecting  a 


l^Inthe  actual  computation* ,  it  l*  advisable  to  utilize  different  color*  to  distinguish  between 
these  values,  the  programmed  flow*  and  the  values  for  the  U.'  vector. 

14  It  would  be  possible  to  consider  these  resistances  as  diffe rent,  depending  on  the  direction 
of  flow.  Uni -directional  flow  could  then  be  assured  by  blocking  out  counter  flow*  with  heavy 
resistance  penalties. 


TABLE  1 

Incidence  Matrix  and  Computing  Routine 
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promising  basis  net  for  the  purpose  of  minimization,  if  thin  set  is  ».ot  optima],  a  series  of 
systematic  iterations  will  secure  such  a  set  in  a  finite  number  of  stages. 

The  basis  thus  initially  rhn«on  js  indicated  by  the  dots  centered  on  the  outside  of  the 
selected  branches  It  consists  of  the  following  branches:  1.  2;  1  6;  2  3;  5.1;  3  4  5,  8;  6,  7; 
and  8,  9;  as  shown  In  Figure  1  1 11  This  set  of  dots  Is  repeated  directly  below  Table  1  in  the  first 
row  labelled  Basis  Indicators.'  It  designates  the  collection  of  column  vectors  which  constitute 
an  active  basis.  Immediately  below  these  dots  the  row  labelled  ’’Program"  shows  the  (unique) 
current  flows  which  are  readily  ascertained  from  Figure  1  by  reference  to  the  branches  in  the 
basis  and  the  influxes  and  effluxes. 

Note  that  57  5  and  2,  5  are  linearly  dependent.  Hence,  at  most,  one  of  these  vectors  may 
enter  Into  a  basis.  The  other  elements  of  a  basis  having  been  chosen,  the  total  amount  of  cur¬ 
rent  to  flow  between  nodes  2  and  5  is  uniquely  determined.  Since  the  objective  is  minimization, 
the  optima]  choice  consists  of  filling  up  each  such  set  of  branches  in  the  order  of  successively 
increasing  Cj's  Thus,  the  flow  of  5  between  nodes  2  and  5  is  allocated  as  follows:  First  as 
much  as  possible  of  this  flow  is  allocated  to  2,  5  Whenever  such  a  link  is  used  to  capacity 
with  overflow  to  one  of  its  parallel  links,  such  as  2.  5,  it  is  not  counted  as  part  of  an  active 
basis  If  For  this  reason  2,  5  receives  a  bar  rather  than  a  dot  in  the  row  of  Table  1  labelled 
"Basis  indicators,"  as  well  as  alongside  the  program  value  4  in  Figure  1.  The  bar  (replacing 
the  dot)  indicates  that  this  branch  is  at  capacity  (with  overflow)  and  is,  therefore,  not  in  the 
basis. 

By  virtue  of  the  linear  dependence  of  the  Kirchhoff  node  conditions,  exactly  eight  vectors 
here  constitute  a  basis.  The  above  method  of  choosing  between  21  and  2,  5  resolves  the  ques¬ 
tion  of  choice  among  links  joining  the  same  node  pair,  only  one  of  which  can  be  in  an  active 
basis 

To  determine  whether  the  programmed  amounts  are  optimal,  the  methods  of  linear 
programming  are  now  applied.  Values,  u>v  are  determined,  one  for  each  node  of  the  network, 
as  follows, 17  One  such  value  may  be  chosen  arbitrarily.  Here  the  first  such  value,  zero,  is 
selected  at  node  1,  and  entered  to  the  left  of  this  node  (under  the  u>  component*')  in  Table  1.*^ 

The  remaining  w  components  are  now  uniquely  determined  by  the  condition  that  each  pair  of 
values  a’j  and  ”,  when  multiplied  by  the  corresponding  elements  of  each  basis  vector  and 
summed,  must  equal  tne  value,  c^,  located  at  the  top  of  each  such  vector.  Call  the  sum  z, 
whira  result  t-om  algebraically  adding  each  pair  of  w's  so  that  the  condition  to  be  satisfied  is 
zj  =  Cj.  Noting  that  zero  appears  as  Wj  opposite  node  1  and  that  cJ2  =  5  it  is  apparent  that 
u’2  -  -5.  Similarly,  since  c16  -  2  it  follows  that  u ^  =  -2  in  order  to  obtain  ui  .w  nlxO 
+  (-1)  (-2)  *  Zj6  =  Cjfi  =  2.  Next,  the  value  u.'2  =  -8  is  determined  from  u>2  -  =  z23  =  c23  =  '*‘ 

Since  u'2  -5  it  is  evident  that  =  -(5+3),  as  shown,  to  achieve  the  desired  equality.  19 

1  5  The  dot  on  branch  5,  4  should  be  ignor,  d  at  this  stage.  It  enters  at  a  subsequent  iteration 
when  3,  4  is  eliminated,  as  will  be  explained  uubsequentl 
i^See  [  3]  or  [7], 

*  7lt  may  be  helpful  to  regard  these  values  as  tentative  node  potentials  in  terms  of  which 
the  flows  are  determined. 

lflNote  that  these  values  may  also  be  determined  directly  from  Figure  1.  Thus,  a  value  of 
zero  is  inserted  at  node  1.  When  c  12  iff  subtracted  from  this  amount,  the  value  of  u>2  at  node  2 

11  secui  ed  Similarly,  the  value  -8  for  node  3  is  obtained  by  subtracting  C11  from  u >7,  and 
so  on.  c *  c 

i^Noti?  that  these  values  may  also  be  determined  directly  from  Figure  1.  Thus,  a  value  of 
zero  Is  inserted  at  node  1.  When  c]2  =5  is  subtracted  from  this  amount  the  value  =  -5  at 
node  2  is  secured.  Similarly  the  value  w3  -8  for  node  3  is  obtained  by  subtracting  c,,  from 
a’j.  and  no  on. 
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Continuing  In  this  fashion  utilizing  only  vectors  designated  by  a  d  >t  bind  therefore  in  the  b.mlsj 
the  r<  m  .lining  vitiu*  s  in  the  left  hand  stub  nf  Table  1  are  thereby  determined  20 

Ail  vectors  in  th’  active  basis,  therefore,  necessarily  have  z.j  -  =  0  If  the  program 

Is  optimal  those  vectors  which  are  not  in  the  active  basis  have  values  ^  -  Cj  02  1  and 
bounded  variables,  such  as  -he  value  for  2,  5.  will  have  values  z^  -  c.  ->0. 

If  anv  value  Is  located  which  does  not  meet  these  conditions,  it  Is  necessary  to  iterate 
further.  In  the  present  instance  z5<1  1  c^4  Tills  vector  is  therefore  brought  into  the  basis. 

as  indicated  by  the  upward  pointing  arrow  just  below  the  first  row  of  "z  values"  of  Table  1 
To  determine  which  vector  Is  to  be  eliminated,  recourse  may  be  had  to  Figure  1  Since  5  4  is 
being  brought  into  the  basis,  the  node  at  4  can  be  secured  via  tills  route.  That  is,  the  vector 

4  may  be  stated  in  terms  of  the  basis  by  tracing  the  unique  path  along  the  basis  tree  either 
from  head  to  tail  or  bail  to  head.  Thus,  5,  4  utilizes  5,  2;  2,  3;  and  3,  4  The  first  is  accorded 
a  minus  sign  as  the  negative  of  2,  5  and  the  latter  are  accorded  positive  signs.  Only  the  latter 
(the  ones  accorded  positive  signs)  are  considered  for  purposes  of  designating  the  vector  to  be 
eliminated  from  the  basis  when  5,  4  is  inserted.  Uy  recourse  to  standard  procedures,  including 
known  methods  nf  resolving  degeneracy.22  the  vector  3.  4  is  eliminated  from  the  basis  and  5,  4 
inserted  instead. 

Crossing  out  tin  program  mud  value  of  /.mi  along  3,  4  and  inserting  it  instead  at  fi,  4, 
all  other  flows  are  left  undisturbed  at  their  previous  values.  Also,  the  only  value  of  which 
requires  alteration,  is  the  one  at  node  4.  Starting  with  w,  =  -7  and  subtracting  c„„  =  3  from 
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this  amount,  the  new  value  u<4  -  -10  is  secured  to  replace  the  old  one  of  -17,  as  shown  in  Fig¬ 
ure  1. 

The  second  row  of  dots  in  Table  1  indicates  the  new  basis  with  the  programmed  amounts 
entered  directly  below  the  dots.  To  check  for  optimality,  all  vectors  not  in  the  active  basis  are 
evaluated  as  before.  The  tost  for  optimality  being  successfully  passed  as  Indicated  by  the 
numbers  In  the  z  value  row  at  the  bottom  of  the  Table.2  '  the  problem  is  solved  with  the  current 
flows  as  Indicated  in  the  "program"  row 

The  methods  for  convex  programming  over  incidence  matrices  here  developed,  may  be 
generalized  to  what  may  be  called  matrices  of  incidence  type  with  nonzero  entries  other  than 
t  1.  spec T'icnUy  they  may  he  considered  as  matrices  for  node  conditions  which  generalize  the 
Kirrhhoff  node  conditions  In  that  the  current  in  a  particular  branch  Is  to  be  multiplied  by  a 
number  associated  with  the  branch  before  entry  into  the  node.  Such  extensions  are  important, 
for  Instance,  In  problems  of  accounting  and  financial  analysis,  and  are  encountered  even  in  very 
simple  cases  such  as  the  ones  described  in  [6],  the  further  refinement  of  these  methods24  to 
the  exact  functions  and  involving  nonlinear  Iteration  s  eps,25  is  also  important  for  dealing  with 
problems  of  risk  and  uncertainty  These  topics  will  be  treated  in  subsequent  papers.25 

2®F:,r  a  discussion  of  the  theory  underlying  these  procedures,  see  [3], 

ZlSince  flows  may  he  ;n  either  direction.it  is  necessary  either  to  use  the  absolute  value  of 
zj  or  else  to  take  account  of  the  direction  of  flow.  In  links  such  as  2,  5,  the  direction  of  flow 
must  be  considered.  The  z.j  which  results  from  a  negative  flow,  is  simply  the  negative  z,  for 
positive  flow.  Thus,  it  is  unnecessary  to  carry  two  columns.  In  either  cas e  —  po sitiv e  or  nega¬ 
tive  flow— on  a  capacitated  link  the  true  z.  must  be  greate r  than  or  equal  to  the  corresponding  C;. 

22See  [2]  J  } 

23l  e.,  |Zj|  <  Cj  for  ill  nonbasis  vectors  and  z,  ,  for  2,  5. 

24bee  footnote  11,  page  1 14.  supra, 

25E  g.,  A.  Charnes  and  fi.  E  Lemkr,  ”'i<  ntinuous  Limit  Methods  in  Mathematical  Pro¬ 
gramming,"  ONR  Research  Mernr.  No.  1,  Systems  Research  Gro  ;p,  Northwestern  University 
( forthcoming) . 
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